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Abstract: We study the Dalitz plot of the Cabibbo-favored charmed-meson decay D~^ 

using dispersion theory. The formalism respects all constraints from analyticity 
and unitarity, and consistently describes final-state interactions between all three decay 
products. We employ pion-pion and pion-kaon phase shifts as input, and fit the pertinent 
subtraction constants to Dalitz plot data by the CLEO and FOCUS collaborations. Phase 
motions of resonant as well as nonresonant amplitudes are discussed, which should provide 
crucial input for future studies of CP violation in similar three-body charm decays. 
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1 Introduction 

Heavy-flavor three-body decays into light mesons provide a valuable source for Standard 
Model tests and beyond. While they are driven, at short distances, by the weak interactions, 
their rich kinematic structure accessible in Dalitz plot distributions makes them a prime 
example for the application of modern tools of amplitude analysis [1]. A major motivation 
for the investigation of heavy-flavor decays is the study of CP violation, which manifests 
itself in the appearance of (weak) phases and requires the interference of different amplitudes 
with, at the same time, different phases in the strong final-state interactions (see, e.g., 
Ref. [2] for an in-depth overview). In contrast to (quasi-)two-body decays occurring at fixed 
total energies, three-body decays offer a resonance-rich environment with rapidly varying 
strong phases throughout the phase space available, which may strongly magnify the effects 
of CP violation in certain parts of the Dalitz plot. 

Obviously, in order to turn the search for potentially very small CP-violating phases in 
such complicated hadronic environments into a precision instrument, it is inevitable to con¬ 
trol the strong dynamics in the final state as accurately as possible, in a model-independent 
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fashion that, however, incorporates a maximum of theoretical and phenomenological con¬ 
straints. The traditional approach to model Dalitz plots in terms of isobars, i.e. a series of 
subsequent two-body decays, and describe the relevant line shapes in terms of Breit-Wigner 
(or Flatte) functions, has clear limitations: it fails to describe in particular the phase mo¬ 
tion of the broad S-wave resonances such as the /o(500) in pion-pion or the iFQ(800) in 
pion-kaon scattering (see e.g. Refs. [3, 4] in the context of heavy-flavor decays), and neglects 
corrections beyond two-body rescattering in an unquantified manner. 

It has therefore been advocated to employ the framework of dispersion theory for ampli¬ 
tude analyses [1], which is built on unitarity, maximal analyticity, and crossing symmetry. 
The dispersive framework adapted to study three-body decays was originally introduced by 
Khuri and Treiman for the decay K ^ 3tt [5] and subsequently further developed [6-10]. 
The formalism has been resurrected in a modern form in Refs. [11, 12]. The Khuri-Treiman 
equations are based on elastic two-body unitarity and explicitly generate crossed-channel 
rescattering between the three final-state particles. The equations are constructed by set¬ 
ting up dispersion relations for the crossed scattering processes, with a subsequent analytic 
continuation back into the decay region. This continuation is performed along the lines of 
the continuation of the perturbative triangle graph and is extensively discussed in Ref. [6]. 

Khuri-Treiman equations have been successfully applied to various low-energy meson 
decays, like e.g. r/ —)• Svr [11-14], a;/(/> —^ Svr [15, 16], or r]' —)• r/Trvr [17]. In this work, we 
extend this formalism to three-body decays of open-charm mesons, analyzing the Cabibbo- 
favored decays D~^ —)• K~TT~^ir^. As input we solely rely on vrvr and ttK phase- 
shift input. While these are not yet decay channels of major interest to study CP violation, 
the final-state interactions are going to be similar for others that are, such as the Cabibbo- 
suppressed decays D —)• 3tt/ttKK. For the decays at hand, inelastic effects are small in large 
regions of the Dalitz plots, and therefore elastic unitarity provides a good approximation: 
the TTvr channel allows for isospin 1 and 2 only, but no isoscalar components, which would 
necessitate a coupled-channel treatment, as a strong coupling to KK occurs. The major 
inelasticities in the ttK channel are found to set in at the r]'K threshold [18-20]. 

Thus with the high-statistics experimental data available [21-23], these decays provide 
a good test case to establish this dispersive framework in higher energy regions and set 
the path to Cabibbo-suppressed decays where traces of physics beyond the Standard Model 
may be searched for. Besides, it allows for a further test of low-energy ttK and tt-tt dynamics 
as well as the importance of crossed-channel rescattering effects in three-body decays. It 
may also provide an insight into scattering phase shifts at higher energies in the future. 

The decay under consideration has been the subject of a number of previous theo¬ 
retical publications, focusing on different issues raised by the experimental results. One 
challenge is the proper treatment of the isospin 1/2 S-wave with the very broad, non-Breit- 
Wigner-shaped A'q(800) (or k) resonance [24], and the inclusion of two scalar resonances 
iFQ(800) and ArQ(1430) in a way that conserves unitarity. Furthermore, the width of the 
iFQ(1430) extracted from the experimental analyses in Refs. [21, 22] is found to be incon¬ 
sistent with PDG values [25]. In addition, the explicit comparison of the ttK partial-wave 
phases extracted from these decays [23, 26] with ttK scattering results [27] seems to indicate 
deviations from Watson’s final-state theorem. 
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Ref. [28] focuses on the isospin 1/2 5-wave final-state interactions, based on coupled- 
channel partial waves for Ktt, Krj, and Krj' constructed dispersively in Ref. [18]. Decay and 
scattering data could be reconciled, although no three-body rescattering effects, isospin 3/2 
components, or ttyt channel were included. Ref. [29] similarly observes mutual consistency 
of ttK scattering and the D-meson decay, using related input to take two-body final-state 
interactions in the ttK isospin 1/2 5- and P-wave into account in terms of the correspond¬ 
ing scalar and vector form factors. Furthermore, the short-distance weak interactions are 
described with the help of an effective Hamiltonian based on a factorization ansatz. Again, 
weak repulsive partial waves (of isospin 3/2 and in the TTTr system) as well as crossed- 
channel rescattering are neglected. We mention that similar approaches, using dispersively 
constructed form factors for two-body rescattering, but neglecting third-particle interac¬ 
tions, have also been applied to P —)• Ktttt decays [30, 31]. 

In Ref. [32], a Faddeev-like equation is solved that builds up three-particle rescattering 
effects. The underlying two-particle ttK amplitudes are obtained form unitarized chiral per¬ 
turbation theory fitted to experimental data. The decay amplitude is simplified to include 
only the isospin 1/2 5-wave, aiming mainly at a study of the importance of rescattering 
effects and the reproduction of the experimental 5-wave phases [23, 26]. The model for the 
weak vertex has subsequently been improved [33]. Ref. [34] applies a similar approach with 
the addition of the isospin 3/2 ttK 5-wave, but is still restricted to 5-waves only. The only 
theoretical analysis known to us with all relevant partial waves, three-particle rescattering 
effects, and effects of the intermediate state included, is Ref. [35]. The author 

performs a full Dalitz plot analysis on pseudo data, which we will later compare to. 

The outline of this article is as follows. Section 2 states some basic kinematical relations 
and shows both isospin and partial-wave decomposition of the decay amplitude in question. 
In Sec. 3, we derive the coupled dispersive integral equations and discuss how to solve these. 
Numerical results are shown in Sec. 4 and compared to experimental Dalitz plot data by 
the CLEO [21] and FOCUS [22] collaborations. We conclude our study in Sec. 5. Some 
technical details are relegated to the appendices. 


2 Kinematics, isospin decomposition, partial-wave expansion 


The Mandelstam variables of the D-meson decay 


D^{pd) A:(pii')7r(pi)7r+(p2) 


( 2 . 1 ) 


are given by s = {pn — Pi)^, t = {pD — P 2 f‘, and u = {pD — Pk)"^■ The corresponding 
scattering angles 6 in the (crossed) scattering processes are given by 


s(t-'u)-A t{s-u)-A t-s 

Zs = COS Us = -, Zt ^ cos Ot = -77;- , Zu ^ cos Ou = 


k{s) ’ ’ ----- 

k{x) = , Ku{u) = 


4M2 

1 -( 2 . 2 ) 


with A = (M|, — ~ ^n) the Kallen function A(x, y, z) = x'^ + y'^ + z'^ — 

2{xy + yz + xz). 
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(a) 


Figure 1. Quark line diagrams of the —>■ 
and internal W~^ conversion (b). 
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(b) 

“7r+7r+/^°7r°7r+ decays: W~^ as a spectator (a) 


We begin with the isospin and partial-wave decompositions of the decay amplitudes 

M _I-+ (D~^ —>• K~Tr~^7r~^) and Adoo+ t X^vr^vr^). We associate the isospin structure 

of the strong hnal-state current in Fig. 1 with the D~^ meson. Since one uu/dd pair is 
strongly produced, the associated isospin of the D meson is given by / = 3/2, 1^ = 3/2. 
Thus the isospin decomposition of the respective (crossed) scattering processes reads 


s / t-channel 

tt-channel 

■b^D+TT°^K0-K+ = ’ 

Md+ko^^Ot,+ = ) 


■Md+K0->-7t+tt0 = ) 

/ , 

-^Z)+iC+—>7r+7r+ ^ ) 


where denotes the amplitude with definite isospin I. The decay amplitudes are given 

by 


■M. —1-+ (s, t. It) T -^D+7r“—>/C“7r+ (T '^) T D+K+^n+TT+ 

Moo+{s,t,u) = Mj:)+^o^KOTr+is,t,u) + Mn+KO^nOn+{s,t,u) + MD+^-^KOno{s,t,u) . 

(2.4) 


We can write down a symmetrized partial-wave expansion simultaneously in s-, t-, and 
ri-channels (the precise relation of which to proper partial waves in a single channel will be 
discussed below). With the expansion in partial-wave amplitudes truncated at the D-wave 
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Figure 2. The associated s-channel scattering diagram via the intermediate states 

. The gray vertex stands for the crossed decay amplitude D+tt” —>■ denoted by 

and the white vertex the —>■ scattering amplitude denoted by i". The dashed line 

gives the contribution to the discontinuity [41], The other channels follow analogously. 


for ttK final states and the P-wave for tt-tt, we obtain 


A4_++(s,t,n) = P^{u) + 


+ [s(t -u)-A] 

+ \ [3(s(^ -u)- Af - k2(s)] (^-^Py^{s) - + (s ^ i)| , 

MQQj^{s,t,u) = Pl{u) + V^{t - s) Fl{u)) + 

+ [s{t -u)-A] ^(s) + ^ [3{s{t -u)-AY- k2(s)] 


- [t{s - n) - A] ( 


- ^ [3(t(s - n) - A)' - ^{t)] , (2-5) 


where the single-variable amplitudes have definite isospin I and angular momentum L in 
the channel associated with the Mandelstam variable featuring as their argument. Note that 
the inclusion of D-waves is somewhat heuristic: in order to rigorously prove the symmetrized 
decomposition (2.5) in the spirit of the so-called reconstruction theorem [12, 36-40], one 
needs to include a subtraction polynomial of higher order (i.e., a larger number of unknown 
parameters) than what we will allow for below. We mainly want to retain the irK D-wave 
to test the effect of the iir|(1430) resonance, which is kinematically accessible in the decay 
phase space. The way we implement this approximately will be discussed in Sec. 3.3. 
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3 Dispersive formalism 


3.1 Unitarity and Omnes solution 


We begin with the dispersive treatment of the associated scattering processes linked to 
the decay by crossing symmetry, D'^it —)• Ktt and D~^K —)• ttyt. The D-meson mass is 
artificially set to < Mk + 2 M 7 r such that the corresponding decay is kinematically 
forbidden. The simpler analytic structure of these scattering processes can be exploited to 
construct dispersion relations for the single-variable amplitudes valid for s,t > {Md + 
and u > {Md + respectively. The analytic continuation back to the physical D- 

meson mass as well as into the kinematic region [Mk + < s,t < {Mjo — 

4M^ < u < {Mjo — yields the anticipated decay amplitudes [6]. 

We demonstrate the framework for the example of the s-channel processes; the t- and u- 
channel amplitudes are constructed analogously. Elastic unitarity gives for the discontinuity 
(see Fig. 2 for A4_|-+) 


2 

discA4_++(s,2;s) = - 
discA4oo+(s,2:s) = ^ 


^ Mij+is, z':)d{l^ - Mf) 5{{q - If - Mf) , 

^ ^ iiJ) 

^ - Mf) 6{{q - If - Mf) , 

^ ’ iij) 

(3.1) 


where T*-^’ ^{x,Zx) —)■ K~7r~^) and T*-^’®°(x, (K*M —>• are the inter- 

mediate-to-final-state scattering amplitudes, q = px + P 2 = ifs, 0) defines the center-of- 
mass frame, in which z^ = cos9f the cosine of the angle between initial and intermediate 
states, and z" = cos9'f the cosine of the angle between intermediate and final state, are 
evaluated. The intermediate-state summation runs over the tuple (i,j) S {(—,-|-), (0,0)}. 
The partial-wave decompositions for the irK (ttyt) amplitudes and full decay ampli¬ 

tudes Mijk read 

I,L 

Mijkis, zf = ^ afff Pl{zs) fl{s) , (3.2) 

I,L 

where the sum runs over isospin and angular momentum components I and L. Further¬ 
more we use the Clebsch-Gordan coefficients o/^i, Legendre polynomials Pl{z), and the 
corresponding partial waves tf{s) and Exploiting the unitarity relation for elastic 

ttK and vrvr scattering we obtain the following partial-wave unitarity relations 

disc/f(s) = 2i fl{s) sin 9{s - sth) , (3.3) 

^Note that in contrast to the definition of the single-variable amplitudes in Eq. (2.5), we have not defined 
the partial waves in Eq. (3.2) to be free of kinematical zeros. This is independent of the singularities these 
partial waves display at the corresponding pseudo-thresholds or upper limits of the physical decay region, 
s = {Md — Mtt)^ or u = {Md — Mk)^, which are well understood, see e.g. Ref. [15] or the discussion in 
Ref. [42] in a perturbative context. 
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where Sj^(s) denotes the elastic hnal-state scattering phase shift. The thresholds in the 
different channels are Sth = hh = {Mk + for ttK and tt^h = 4:M^ for tttt scattering, 

respectively. Since the discontinuity of f[ and the according single-variable amplitude 
coincide on the right-hand cut, we have 

disc fl{s) = disc J-'f (s) 

^ fLis) = {His) + His)) , ( 3 . 4 ) 

where we have introduced the inhomogeneities P^is) that are free of discontinuities on the 
right-hand cut by construction. They incorporate the left-hand cut contributions and will 
be further discussed in Sec. 3.2. From Eqs. (3.4) and (3.3) we obtain 

discT'f (s) = 2i {His) + Piis)) e{s - sth) sin , (3.5) 


which has the form of an inhomogeneous Hilbert-type equation. The homogeneous solution 
P[is) = 0 is given by the so-called Omnes function H£(s) [43] times an analytic function 


Pfi 


His) = p[{s)ni{s) , ni{s) = expj- f 

Pth 

The inhomogeneous solution is obtained by a product ansatz 


s'(s'-s) 


(3.6) 


His) 




ds' sin(5£(s').T^(s') 1 
s'” I (s') I (s' - s) / ’ 


(3.7) 


where P[is) is now a polynomial of order n — 1, and the number of subtractions n is chosen 
such that the convergence of the dispersion integral is guaranteed. 

As our approach relies on elastic unitarity (see Ref. [44] for a generalization of the 
Khuri-Treiman formalism to coupled channels), the formalism breaks down when inelastic 
channels become important. We assume that Watson’s theorem [45] is a good approximation 
up to the r]'K threshold in the ttK channel. Inelastic effects in the prominent ttK S'-wave 
systems are found to become sizable above the rj'K threshold [18-20]. The main inelastic 
contributions in the isospin 1/2 P-wave come from the ttK* and pK channels, which become 
noticeable in the energy region where they couple to AT* (1410) and AT* (1690) [20]. In all 
exotic partial waves, i.e. the isospin 2 vrvr system as well as / = 3/2 vrAT partial waves, 
inelastic effects are assumed to be negligible. 


3.2 Inhomogeneities 


With the scattering phase shifts given as fixed input, the only quantities left in the dispersion 
integrals Eqs. (3.7) are the inhomogeneities which are determined as the projections of 
the crossed-channel amplitudes onto the considered channel. They can be re-expressed in 
terms of the single-variable amplitudes pj^{x), such that we obtain integral equations that 
can be solved for the p[(x). With the aid of Eq. (3.4) we find 


flip = 


2L + 1 


2a 


ijk 


dzxMl^f^ix, Zx)PLizx) = Hi^ix){Pl{x) +P[{x)) 


I,L 

2L + 1 


= ijk ^ / dZxMljf.ix, Zx)PLiZx) - Plix), 

ZQjJ Kj yX j J —1 


(3.8) 


7 






where A4jjf,(x, Zx) denotes the projection of the full decay amplitude Mijk[x,Zx) onto 
isospin Ix eigenfunctions in the x-channel. One term of the projection integral over 
^x) is always F^{x), such that the right-hand-cut discontinuity is canceled. The 
inhomogeneities are then indeed free of discontinuities on the right-hand cut as anticipated. 
The resulting inhomogeneities are given in Appendix A. 

The interpretation of Eq. (3.8) as an angular integration is valid in the scattering region 
and needs to be analytically continued into the unphysical and decay regions. Performing 
the angular integration naively in the decay region results in crossing the unitarity cut. 
The prescription on how to perform the continuation has been extensively discussed in 
Ref. [6], motivated by the continuation of the (perturbative) triangle graph into the decay 
region. It ultimately leads to the prescription —)• M|, -|- ie, which allows one to derive 

an integration path that avoids the unitarity cut. 


3.3 Number of subtraction constants 

The minimal number of subtractions needed is dictated by the asymptotic behavior of the 
integrands in Eqs. (3.7). The decay amplitude and thus the inhomogeneities are assumed to 
grow at most linearly asymptotically, loosely based on the Froissart bound [46]. Assuming 
the phase shifts to approach constant values Sj^(oo) for large energies, the Omnes functions 
fl£(x) behave like oc asymptotically. With the following assumption for the phase 

shifts 5^ at high energies: 

lim 5y^(x) = 2tt , lim 6V‘^{x) = vr, lim 5y'^{x) = vr, 

X^OO X^OQ 


lim (x) = 0 , 

C^OO ^ 

lim (5 q (x) = 0 , 
x—>-oo 


lim (5y^(x) = 0 , 
c—>-oo 

lim 5}(x) = vr, 

ir^oo 


lim 6^^ 

X^OO 


(x) = 0, 


(3.9) 


we need two subtractions for J^q, and four subtractions for and one sub- 

traction for to obtain convergent dispersion integrals. Note that the difference in the 
number of subtractions for and despite identical phase asymptotics, is due to the 

different kinematic prefactors for P-waves with equal and unequal masses, see Eq. (2.5). 

needs no subtraction, but as the ttK isospin 3/2 P-wave phase shift is very small and 
assumed to vanish at high energies, we neglect it altogether. Similarly, also the I = 3/2 
D-wave is put to zero. 

The inclusion of the D-wave T 2 is delicate. Formally it requires no subtractions, 
but the kinematical pre-function corresponding to the L = 2 Legendre polynomial, multi¬ 
plied with the required momentum factors to make it free of kinematical singularities, see 
Eq. (2.5), violates the assumed high-energy behavior of the decay amplitude and thus of 
all inhomogeneities. Therefore we will follow a “hybrid approach” for the D-wave: we will 
only consider the projections of S- and P-waves of other channels in order to generate the 
D-wave inhomogeneity, but will exclude D-wave projections, thus eschewing the need for 
further subtractions. This is loosely motivated by analogous observations in low-energy pro¬ 
cesses calculated in chiral perturbation theory, where higher partial waves are dominated 


by crossed-channel loop diagrams that correspond to low partial waves in those crossed 
channels. 

In total we have eleven subtraction constants. However, the resulting representations 
of the decay amplitudes Eqs. (2.5) are not unique due to the linear dependence of the 
Mandelstam variables s, t, and u: One can construct polynomial contributions to the 

single-variable amplitudes that leave the complete decay amplitudes A4 _and 

A^00-i-('®’invariant; this is obvious in a standard dispersive representation, however 
slightly less trivial to demonstrate in the Omnes representations discussed above [47]. The 
polynomial coefficients can be tuned such that a maximal number of subtraction constants 
is eliminated to obtain a linearly independent set. These polynomials span the so-called 
invariance group of the decay amplitudes. Details are discussed in Appendix B. We choose 
to eliminate the subtraction constants in the nonresonant 7 = 3/2 ttK and 7 = 2 Trvr S- 
waves, the rationale being solely to retain them in presumably large, resonant partial waves. 
This leaves seven linearly independent complex subtraction constants, 


T^{u) 

Tl{u) 




du' {u') sin 6q {u') 


TT 


'«th - u) ’ 

. f V? du' P}{u')s\n5\iu') 

co-7ciuH-/ —= I ^ 




q'^^(s)<[ C2 -I- C3S + C4P + C5P + 


n: 


f(«) 


^ 7°° dg^ /o^^(s0sin,5g^^(s0 

^ Ah 

ds' (s') sin 5 }^^ (s') 
|d}'^^(s')|(s'— s) 


s 

TT 

3/2/ 


poo 


-u) 

ds' pQ^'^(s') sin (5^^ (s') 


e/4 




1/2, 


s')|(s'- 


C6 


poo 


D^/'(s) 


,7^2''^(s') sin 52^^(s') 


1/2, 


ds' 


TT 


Sth 


\n 


1/2, 


3 ')|(S' - 


(3.10) 


The subtraction constants cannot be determined in the framework of dispersion theory and 
have to be obtained either by matching to a more fundamental dynamical theory, or, as in 
this work, by a fit to experimental data. The solution space of the coupled system Eq. (3.10) 
has thus dimension seven, corresponding to the seven complex subtraction constants. Since 
the equations depend linearly on the subtraction constants, it is convenient to choose seven 
independent basis sets and solve the equations for each of these sets. We call those solutions 
basis functions. In particular, we choose for the ith basis function Aii(s,t,u) the set of 
subtraction constants Cj = 6ij with i,j = 0... 6. The full solution Ai(s,t,u) is then 
obtained by 

M.(s,t,u) =''^^CiM.i(s,t,u). (3-11) 

i 

The basis functions are entirely determined by the phase shift input, as well as the masses 
of all particles involved (taken from Ref. [25]). The phase shifts are obtained from solutions 
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of the vrvr Roy equations by both the Bern [48-50] and the Madrid-Krakow [51] groups, as 
well as the Roy-Steiner equations for ttK scattering solved by the Orsay group [52]. 

3.4 Solution strategy 

In this section we discuss different solution strategies of the Khuri-Treiman-type equa¬ 
tions (3.10), their issues, and present our new solution strategy. 

The standard solution strategy for the linear coupled double integral equations (3.10) 
has been an iteration procedure as performed for example in Refs. [13, 15, 17] or numerically 
faster with the introduction of integral kernels in Ref. [53]. Starting with an arbitrary in¬ 
put for the single-variable amplitudes (e.g. just the Omnes functions), the inhomogeneities 
are evaluated; with these the dispersion integrals are determined to obtain a new set of 
single-variable amplitudes. This cycle is repeated until satisfactory convergence is reached. 
Unfortunately, the convergence of this iterative procedure is not always guaranteed, de¬ 
pending on the mass of the decaying particle and the number of subtractions: for larger 
decay masses and more subtractions applied, the corrections in each iteration step can be 
too large to reach the hxed-point solution. We hnd this to be the case in the H-meson 
decays considered here. 

This necessitates a different solution strategy. Since the set of integral equations is 
linear in the single-variable amplitudes it is convenient to set it up in the form of a matrix 
equation instead. Provided that the matrix is invertible a unique solution exists. One 
such inversion strategy is known as the Pasquier inversion [9] (see Ref. [54] for a recent 
comparison of Pasquier inversion and iterative solution), where a method to reduce the 
double integral equation to a single integral equation is introduced. The procedure involves 
the deformation of the integral contours of both integrals, allowing one to interchange the 
order of integrations such that a unique kernel function is obtained. The coupled single 
integral equations thus obtained do allow for a direct solution via matrix inversion. 

We will follow a slightly modified strategy, constructing a matrix equation without per¬ 
forming a Pasquier inversion. In this context it is beneficial to solve for the inhomogeneities 
instead of the single-variable amplitudes, the advantage being that the inhomogeneities need 
to be evaluated only on the right-hand-cuts. The single-variable amplitudes themselves can 
be obtained in the whole complex plain in a straightforward manner by performing the 
dispersion integral over the inhomogeneities once. 

To illustrate the solution strategy we limit ourselves to one hypothetical inhomogeneity 
equation without any loss of generality, 



(3.12) 


and focus on the functions J-{s) = that are free of singularities at the 


pseudo-threshold or upper limit of the kinematically accessible decay region (which is a 
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zero in k(s)). Inserting Eq. (3.7) into Eq. (3.12) yields 


/(.) = 


.2L+1 


(^) 


+ 


j d2:sz™0(t(s,z^))|p(t(s,z^)) 

'■ r°° dx T{x)sm6{x) 

J Si 


2 

t{s,zsT_ rdx 

*Sth 


x"- |r2(x)|K2^+l(x)(x — t(s, 2;*)) 

1 ~ 

= A(s) H— / F{x)K{s,x)dx . 

^ J Sth 


(3.13) 


The function ^(s) contains the dependence on the subtraction polynomial, while the inte¬ 
gration kernel K{s,x) is independent of any subtraction constants: 


7l(s) = 


^2L+lfg\ 


K{s, x) = 


'-1 


dzs z^ p{t{s, Zs)) Q{t{s, Zs)) , 


sin S(x) 


x'^\Q-{x)\k^^+^{ x) J_i 


\ t{s,ZsT zfn{t{s,Zs)) 


X - t{s,Zs) 


(3.14) 


Equation (3.13) is thus a linear integral equation for F{s), to be solved for a given set of 
subtraction constants. Discretizing Eq. (3.13) yields 


Msi) =(^^ij - K{si,x)dx'^P{sj), 


(3.15) 


which is solved by matrix inversion; the numerical treatment is relegated to Appendix C. 


4 Numerical results and experimental comparison 

Solving the coupled integral Eqs. (3.10) with the algorithm presented in the previous section, 
we obtain the single-variable basis functions (.T^)i depicted in Figs. 3 and 4. The vector 
resonances 77* (892) (in and p(770) (in pj) as well as the ttK D-wave resonance 

i7|(1430) (in are clearly visible. The basis functions include the effects of the 

scalar states i7Q(800) and i7g(1430), while the exotic Tq and basis functions are free 
of resonances. 

The error bands in Figs. 3 and 4 are determined by a conservative error estimate of the 
phase shifts: For the 5-wave 'kK and tttt phases the error is assumed to rise linearly from 
zero at threshold to ±20° at 2 GeV. Beyond 2 GeV the error is fixed to ±20°. The ttA 
isospin 1 /2 P- and D-wave phase errors and tttt P-wave phase errors are similarly obtained, 
with the only difference that the linear rise of the error sets in after the 77* (892), 772(1430), 
and p{770) resonances, respectively. In the ttyt P-wave case we additionally vary between 
the phase-shift data from Refs. [48-51]. 

In the following we compare our theoretical decay amplitude to the experimental D~^ —)• 
77 “7r+7r^ Dalitz plot data from the GLEO [21] and FOCUS [22] collaborations. Exploiting 
the symmetry of the process under the interchange of the two pions, we can restrict the 
comparison to the region s < t by mirroring the remaining half of the Dalitz plot into this 
region. 
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[GeV] 


Figure 3. Real (red) and imaginary (blue) parts of the single-variable functions for i = 

0, ..., 3. The vertical dashed lines denote the kinematical limits of the decay region. 


The experimental events are collected in equidistant bins of size 0.044 GeV^ x 
0.044GeV^. Bins which overlap with the phase space boundary are discarded, resulting 
in 493 bins over the considered fit region {s < t < {Mk + M^iY). The following event 
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Figure 4. Real (red) and imaginary (blue) parts of the single-variable functions for i = 

4, ..., 6. The vertical dashed lines denote the kinematical limits of the decay region. 


distribution function was used for the fit analogously to the experimental analyses 




rSi-\-S - 

J ti — S J S-i — S 


UgAfs\M-++{s,t,u)fe{s,t) + {l-Ug)AfBB{s,t) dsdt, (4.1) 


with being the center of the corresponding bin and 26 the bin width, e(s,t) the 
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efficiency parametrization, B{s,t) the background parametrization, A4ig and Mb normal¬ 
ization constants such that the background and signal term are normalized to unity, and 
the signal fraction /sig. 

We minimize the following 


2 ^ - (#events/bin)i 

X =2^ — 


i=0 


(^events/bin)j 


(4.2) 


where M is the number of events, the sum runs over the number of bins and the error 
on the binned data is assumed to be purely statistical. In addition to the full dispersive 
representation Eq. (3.10), we also ht a simplihed decay amplitude to data, which is given 
by a sum of Omnes functions multiplied by polynomials: 


+ [■s(i - u) - A] Sd2^(s) ^ 3{s{t -u)-AY- k^{s) 

-|- (t s), 


a 


M) 

(4,3) 


where the c( are again complex fit constants. Equation (4.3) emulates a dispersively im¬ 
proved isobar model that neglects any crossed-channel rescattering effects. The number of 
polynomial fit constants is chosen to resemble the number of degrees of freedom in the full 
dispersive result Eq. (3.10) as far as possible; with certain caveats that preclude an imme¬ 
diate quantihcation of three-particle rescattering effects in the same straightforward way as 
performed for dvr decays in Ref. [15]. In Eq. (3.10), two subtraction constants cq and ci 
are contained in the vrvr P-wave, which only contributes indirectly via the intermediate state 
K^7r^7r~^ to the decay and thus does not show up in the pure Omnes amplitude Eq. (4.3). 
In addition, every Omnes function in Eq. (4.3) needs at least a normalization constant to 
adjust the strength of individual amplitudes, while some single-variable amplitudes do not 
have any subtraction constants. Einally, once the D-wave is included we have one additional 
complex fit parameter Cy in the pure Omnes fits. For that reason we consider both Omnes 
and the full dispersive hts without (Omnes 1, full 1) and with D-wave (Omnes 2, full 2). 

We have the freedom to fix one subtraction constant, as both the overall normalization 
and the overall phase are arbitrary and factorized out; we choose C 2 = C 2 = 1. This leaves 
13 (15) real fit constants for the full / Omnes fits. 

Following experimental custom, we will employ so-called fit fractions to characterize the 
relative importance of various single-variable functions. These are defined in the following 
way 

pp/ ^ J {Pjjxjs, t)) Pjjxjs, t))p ds dt 
J\Ml _|_+(s, t, ri)p ds dt ’ 

where the Pj denote the angular prefactors of the corresponding single-variable amplitudes 
in the total amplitude. The integration runs over the fitted Dalitz plot region. In general 
these fit fractions are not unique due to the freedom of adding an element of the invariance 
group Eq. (B.l); the projections onto partial-wave amplitudes then will lead to different fit 
fractions. 
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Full 1 

Full 2 


Omnes 1 

Omnes 2 

Co X GeV^ 

2.7 ±0.8 

1.2 ± 0.2 

Ic'ol 

0.9 ±0.3 

0.9 ±0.7 

|ci| X GeV^ 

3.8 ± 1.2 

2.2 ±0.5 

|c'il 

3.0 ± 1.5 

4.0 ± 1.3 

C 2 

1 (fixed) 

1 (fixed) 


1 (fixed) 

1 (fixed) 

C3 X GeV^ 

2.8 ±0.4 

2.2 ± 0.1 

Ic^jl X GeV^ 

1.9 ±0.2 

2.0 ± 0.2 

|c 4 | X GeV^ 

2.0 ±0.5 

1.4 ±0.1 

|c^| X GeV^ 

0.9 ±0.1 

1.1 ± 0.1 

[csl X GeV® 

0.7 ±0.3 

0.4 ±0.1 

|c(i| X GeV® 

0.13 ±0.3 

0.19 ±0.02 

[cel X lO^GeV^ 

4±3 

6±2 

|c(i| X GeV^ 

0.11 ±0.05 

0.10 ±0.03 




lc '71 X lO^GeV® 

— 

6±4 

argco 

0.1 ± 0.2 

1.1 ±0.3 

argc() 

0.2 ± 0.8 

0.4 ±0.4 

arg Cl 

0.3 ±0.2 

1.2 ±0.3 

arg c'l 

-0.8 ±0.3 

-0.4 ±0.2 

argcs 

- 0.2 ± 0.1 

0.0 ± 0.1 

arg C 3 

0.2 ± 0.2 

0.3 ±0.2 

argC4 

-0.5 ±0.1 

0.0 ± 0.1 

argC4 

0.4 ±0.2 

0.2 ± 0.2 

argcs 

- 0.1 ± 0.1 

0.1 ± 0.1 

arg C 5 

0.2 ±0.4 

0.0 ±0.3 

argcg 

-0.3 ± 1.2 

-0.9 ±0.2 

arg Cg 

0.0 ± 0.1 

0.0 ±0.3 




arg c'y 

— 

0.4 ±0.3 

xVd-o.f. 

1.18 ±0.03 

1.10 ± 0.02 


1.30 ±0.06 

1.08 ± 0.02 


Table 1. Fit to CLEO data: Numerical fit results for the subtraction constants Ci and c' and the 
corresponding x^/d.o.f.. Four fit scenarios are considered: the full dispersive fit, without D-wave 
(full 1) and with D-wave (full 2), and the Omnes fits of Eq. (4.3), without D-wave (Omnes 1) and 
with D-wave (Omnes 2). The errors on the parameters are evaluated by varying the basis functions 
within their error bands. 


4.1 Comparison to the CLEO data 

The Dalitz plot measured by the CLEO collaboration [21] contains 140793 events. The 
efficiency and background parametrizations are given explicitly.^ Our fit results are sum¬ 
marized in Table 1, together with the fit fractions in Table 2. In the full dispersive fits (full 
fits 1/2), the resulting values for the subtraction constants in Table 1 have similar order 
of magnitude with the exception cg, which is rather small. This can be understood by the 
large single-variable amplitude in this particular basis function (see Fig. 4). Further- 

^ The threshold factors T{x) used in there read [55] 

rp( ) — I {'n'Eth,x\x 3^max I), for 0 < Eth,x\x ^max I < 1/2, 

I 1, for Eth,x\x - XmaxI > 1/2 . 
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Fit 

FF 2 

2 X FF^^ 

2 X 

2 X FFq/^ 

2 X FF 2 ^^ 

Full 1 

(37 ±23)% 

(190 ± 60)% 

(11 ± 3)% 

(65 ± 35)% 

— 

Full 2 

(8 ± 3)% 

(72 ± 12)% 

( 10 ± 2 )% 

(16 ± 3)% 

(0.1 ±0.05)% 

Omnes 1 

(48 ± 16)% 

(178 ± 22)% 

(7 ± 1)% 

(395 ± 35)% 

— 

Omnes 2 

(9.5 ± 8 )% 

(91 ± 22)% 

(8 ± 0.5)% 

(240 ± 40)% 

(0.13 ±0.03)% 


Table 2. Fit fractions CLEO: The resulting fit fractions of Eq. (4.4) for the different fit scenarios; 
the errors on the parameters are evaluated by varying the basis functions within their error bands. 
The fit fractions for the ttK amplitudes are multiplied by two to account for the s t symmetry. 


more the phases of the subtraction constants (cq, ci) nearly agree modulo vr. The same 
holds for the Fq subtraction constants (c 2 to C 5 ) especially for the full fit 2. This suggests 
that with overall phases factorized, the subtraction constants for the F\ and likewise the 
Fq amplitude are almost real. The differences of the single-variable amplitude phases 
to the elastic phase shifts depicted in Fig. 7 are thus predominantly due to the dispersion 
integrals, i.e. the crossed-channel rescattering effects. 

Including the D-wave improves the y^/d.o.f. slightly from 1.18 ± 0.03 to 1.10 ± 0.02. 
Note that in the full dispersive representation, no additional fit constants are introduced 
when the D-wave is added. The inclusion of the D-wave does not change the phases of most 
subtraction constants beyond their uncertainties, with the exception of cg; the magnitudes, 
in contrast, change signihcantly for almost all subtractions. Considering the fit fractions 
in Table 2, we observe that the inclusion of the D-wave in the full ht 2 reduces the highly 
destructive interference between the two S-wave amplitudes in the ttK channel. We wish 
to point out that also in Ref. [22], a large cancellation between the isospin 1/2 and isospin 
3/2 S-wave components of —164% is seen, with individual ht fractions of (207 ± 24)% and 
(40 ± 9)%, respectively, which show a comparable behavior to our full ht 1. Although the 
ht fraction of the D-wave itself is very small, it thus has a rather large impact on the S- 
and P-waves. A similar phenomenon is seen in Ref. [21] where the ht quality deteriorates 
considerably when removing the small D-wave. Although we do not ht the whole Dalitz plot, 
the ht fractions for the resonant single-variable amplitudes for fI^"^ and F^^'^ agree 

well with the results from Refs. [21, 22]. The Fq ht fraction corresponds to the isospin 2 
TTTT S-wave component of FF ~ (9.8... 15.5)% found in Ref. [21] within different ht models, 
and together with the ht fraction of F^^"^ agrees with the nonresonant contribution found 
in Ref. [22] of FF (29.7 ± 4.5)%. 

Although the Omnes hts (Omnes 1, 2) yield overall similar results, the strengths of 
the individual amplitudes shown in Table 2 are highly implausible and probably sufficient to 
reject this model. In particular the contribution of the nonresonant isospin 3/2 ttK S-wave 
is vastly beyond all reasonable expectations, and cannot be justihed. In contrast to the full 
ht, this situation is not ameliorated signihcantly by including the D-wave. We conclude 
that crossed-channel rescattering effects are essential to obtain sensible ht fractions. 
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s [GeV^] 



Figure 5. From top to bottom: The experimental data from CLEO [21] depicted in a binned Dalitz 
plot. Below that the theoretical Dalitz plot fitted to the data (fit 2). The dashed line denotes 
the rj'K threshold. The lowest plots show slices through the Dalitz plot. The red and blue curves 
correspond to the full fits 1 and 2, respectively. 
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Fit constant 

Full 1 

Full 2 


Omnes 1 

Omnes 2 

Co X GeV^ 

3.0 ±0.8 

0.6 ±0.3 

ko 

0.4 ±0.2 

0.6 ±0.3 

|ci| X GeV^ 

3± 1 

0.9 ±0.3 

Ic'il 

1.9 ±0.8 

2.2 ±0.5 

C2 

1 (fixed) 

1 (fixed) 

4 

1 (fixed) 

1 (fixed) 

|c3 X GeV^ 

2.8 ±0.8 

1.9±0.1 

Cg X GeV^ 

1.7 ±0.2 

1.8 ±0.2 

|c4| X GeV^ 

2.5 ±0.6 

1.1 ±0.1 

|c^| X GeV^ 

0.9 ±0.2 

1.0 ±0.2 

Icsl X GeV® 

0.4 ±0.2 

0.3 ±0.1 

Ic'sl X GeV® 

0.1 ±0.2 

0.3 ±0.1 

Icel X GeV^ 

0.2 ±0.1 

0.0 ±0.1 

Ic'el X GeV^ 

0.1 ±0.4 

0.1 ±0.1 




jcf,] X lO^GeV® 

— 

7±4 

argco 

0.5 ±0.3 

0.9 ±0.3 

argCo 

0.7 ±0.5 

0 ± 1 

arg Cl 

0.6 ±0.4 

1.1 ±0.2 

arg c'l 

-1.1 ±0.4 

0.2 ±0.3 

argcs 

0.0 ±0.2 

0.0 ±0.1 

argcg 

0.4 ±0.2 

0.2 ±0.2 

argC4 

-0.2 ±0.3 

0.0 ±0.1 

argC4 

0.6 ±0.2 

0.2 ±0.3 

argcs 

0.2 ±0.3 

0.0 ±0.1 

arg Cg 

0.8 ±0.2 

0.2 ±0.3 

argce 

-0.6 ±0.7 

-1.0 ±0.4 

argcg 

-0.7 ±0.3 

-0.9 ±0.3 




arg Cj 

— 

-1.1 ±0.5 

xVd-o.f. 

1.20 ±0.01 

1.21 ±0.02 


1.25 ±0.02 

1.17 ±0.01 


Table 3. Fit to FOCUS data: Numerical fit results for the subtraction constants Ci and c' and the 
corresponding x^/d.o.f.. The same four fit scenarios as in Table 1 are considered. The errors on 
the parameters are evaluated by varying the basis functions within their error bands. 

The resulting Dalitz plot as well as a one-dimensional representation in terms of slices 
through it are displayed in Fig. 5. The bin numbering for the latter is organized in terms 
of t-slices for constant s, subsequently glued together with the next slice of higher s. We 
evaluate the event distribution function Eq. (4.1) over each bin and compare to experimental 
data. The rather small error band on the fit results suggests that the uncertainty in the basis 
functions is largely compensated by interference effects between the different single-variable 
amplitudes, as well as by corresponding variations in the fitted subtraction constants. 

4.2 Comparison to the FOCUS data 

The FOCUS Dalitz plot data [22] includes 52460 ± 245 signal and 1897 ± 39 background 
events. With the resulting signal fraction of ~ 96.5% we perform the full and Omnes fits 
as above. Table 3 summarizes the fit results together with the fit fractions in Table 4. The 
overall picture is very similar to the CLEO fit results with a slightly bigger y^/d.o.f. ~ 1.2. 
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Fit 

oto 

2 X FF^^ 

2 X FfJ/^ 

2 X FFq/^ 

2 X FF2^^ 

Full 1 

(12±4)% 

(59 ± 25)% 

(7.5 ± 2.5)% 

(39 ± 27)% 

— 

Full 2 

(5 ± 3)% 

(67± 10)% 

(12±1)% 

(8 ± 6)% 

(0.17 ±0.07)% 

Omnes 1 

(33 ± 17)% 

(91 ± 37)% 

(9 ± 1)% 

(215 ± 135)% 

— 

Omnes 2 

(89 ± 42)% 

(20 ± 12)% 

(11 ± 1)% 

(180 ± 60)% 

(0.4 ±0.05)% 


Table 4. Fit fractions FOCUS: The resulting fit fractions of Eq. (4.4) for the different fit scenarios; 
the errors on the parameters are evaluated by varying the basis functions within their error bands. 
The fit fractions for the ttK amplitudes are multiplied by two to account for the s t symmetry. 


The Omnes fits again result in nonphysical fit fractions (see Table 4), and from here on we 
will only compare the full fits of both experimental data sets. Starting with the fit without 
D-wave (full fit 1) we observe similar moduli of the subtraction constants compared to the 
CLEO results, however the phases do differ. The fit does not show the large destructive 
interference effects between the isospin 1/2 and isospin 3/2 S-wave that we find in the 
CLEO ht. 

No improvement in the y^/d.o.f. is observed when we include the D-wave (full fit 2). 
However the contribution from the nonresonant amplitudes, the isospin 2 and isospin 3/2 
S-waves, are reduced (see Table 4). The fit fractions of the full fit 2 differ slightly from the 
CLEO hts; in particular the nonresonant 5-waves contribute less in the FOCUS data. 

In the full fit 2 the phases of the Jq subtraction constants persist to nearly agree 
modulo vr; the same holds for T/' subtraction constants. It is reassuring that the overall 
picture of the phases of various subtraction constants is consistent in the full fit 2 results 
for both CLEO and FOCUS. 

In Fig. 7, we compare moduli and phases of the resulting single-variable amplitudes 
as htted to the two data sets; the phases are also compared to the input phase shifts used 
in the Omnes functions. The resulting phase motions largely agree in the two analyses 
within uncertainties, with the possible exception of some deviations in Dq in the region 
of the Kq (800) resonance, where the phase extracted from the CLEO ht rises more quickly. 
There are signihcant deviations from the input phase shifts throughout: there is no naive 
realization of Watson’s theorem in the presence of three-body rescattering effects, see e.g. 
recent discussions in Refs. [54, 56]. This is also the explanation for the observed discrepancy 
of the ttK 1 = 1/2 S-wave phase as extracted from these decays by the E791 [23] and 
FOCUS [26] collaborations, compared to the scattering phase-shift analyses [27]: while the 
phase shift rises to about 67°-97° at y/s = 1.3 GeV [52], the experimental analyses of D- 
decay data suggest an increase in the phase from threshold by about 133°-164° (read off 
via Ref. [32]). Figure 7 shows that in the dispersive formalism, the phase at 1.3 GeV is 
about 182°-198° (CLEO) or 170°-183° (FOCUS)—even larger than found in Refs. [23, 26]. 
We emphasize that these results are based on a formalism that uses the scattering phase 
shifts [52] as input: the deviations in the decay amplitude S-wave are due to complex phases 
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s [GeV^] 



s [GeV^] 



bin number 


Figure 6 . From top to bottom: The experimental data from FOCUS [22] depicted in a binned 
Dalitz plot. Below that the theoretical Dalitz plot htted to the data (fit 2). The dashed line 
denotes the rj'K threshold. The lowest plots show slices through the Dalitz plot. The red and blue 
curves correspond to the full fits 1 and 2, respectively. 
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Figure 7. Left column: Absolute values of the single-variable amplitude in arbitrary units of full 
fit 2 (CLEO in red, FOCUS in blue). The overall normalization is chosen such that the absolute 
values in the AT* (892) peak agree. Right column: Phases of the single-variable amplitudes (CLEO: 
red, FOCUS: blue) and input scattering phases (black) in radiant. The phases are fixed to zero 
at the two-particle (tttt, ttK) thresholds. The dotted lines visualize the fitted area; for the ttK 
amplitudes from threshold to the rj'K threshold and the full phase space for the tttt amplitudes. 


induced by three-body rescattering effects. 

In general, the corrections compared to input phase shifts are smallest for narrow 
resonances, in particular in the 7=1/2 ttK P- and D-waves. The largest phase differences 
are observed in the nonresonant amplitudes, where the phases of Pq and show a 27r rise 
due to zeros in imaginary or real parts close to threshold in individual basis functions. Note 
how these seemingly drastic differences are accompanied by very small absolute magnitudes 
of the amplitudes in question: in view of the aim to control the phase behavior of the 
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complete, combined decay amplitude accurately, these specific deviations are still rather 
small. 

Turning to the moduli of the single-variable amplitudes, the relative strength of the 

K*{892) resonance (in compared with the K|(1430) (in agrees between CLEO 

1 /2 

and FOCUS hts. However the dip in the J-'q' amplitude is shifted to higher energies in 
the FOCUS ht and slightly more pronounced. The moduli of the nonresonant amplitudes 
and Tq turn out to be smaller in the FOCUS fit, which is also underlined by the fit 
fractions (compare Tables 2 and 4). 


4.3 Comparison to other approaches 

So far the only theoretical approach known to us that includes all relevant partial waves, 
three-particle rescattering effects, and the isospin coupled intermediate state is 

Ref. [35]. The treatment is based on a unitary coupled-channel framework. The two-particle 
rescattering contributions are hxed by the ttRT and ttyt scattering data, phases, and moduli. 
Three-body rescattering effects are generated by solving a Faddeev equation. In addition 
to the three-body rescattering a three-body potential, based on hidden local symmetry, is 
introduced modeling vector meson exchanges. The author studies the influence of individual 
rescattering contributions by considering different fit scenarios; crossed-channel rescattering 
effects and three-body potential turned off (isobar fit), three-body potential turned off (Z 
fit), and the full fit. An additional contact term breaking unitarity is allowed for, which in 
the full ht turns out to be negligible. The decay amplitude depends on 27 to 39 degrees 
of freedom depending on the considered ht model, which is more than twice the number of 
parameters included in our full ht. 

To compare the ht fractions obtained in Ref. [35], we note that the isobar ht theoret¬ 
ically compares closest to our Omnes hts, while the Z ht does to our full hts. However 
the isobar ht has a large contribution from the unitarity-breaking contact term (considered 
as a “background” contribution) of 17.7%, such that a direct comparison is not sensible. 
Concerning the full and the Z ht, a large destructive interference between the isospin 1/2 
and isospin 3/2 S-waves is seen, similar to our CLEO ht 1 conhguration. The isospin 1/2 
P-waves are of similar size, ~ 15% compared to our 10 — 14%, but the Trvr P-wave contri¬ 
bution is smaller (1.8 — 3.8%) than our contributions in either full ht 1 or CLEO full ht 2. 
It agrees only with the FOCUS full ht 2. — Concerning this comparison, we should stress 
once more that in contrast to Ref. [35], we do not ht the full Dalitz plot. 

Unfortunately the improvement due to crossed-channel rescattering cannot be quanti- 
hed in a simple way in Ref. [35] either. The improvement going from the isobar to the Z 
and then further to the full model can also be due to the introduction of further degrees 
of freedom; as discussed above, we encounter a similar problem in our analysis. However 
the background term, which gives an indication for missing physics, reduces dramatically 
once the crossed-channel rescattering effects and the coupled intermediate state 
are included. This is a similar conclusion as drawn from the dispersive analysis oi ^ ^ Zn 
Dalitz plots [15], which rendered phenomenological contact terms [57, 58] superhuous. 
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5 Conclusion 


In this paper we have analyzed the —>• decay with a dispersive framework 

based on the Khuri-Treiman formalism that satisfies analyticity, unitarity, crossing sym¬ 
metry, and includes crossed-channel rescattering among the three final-state particles. 

The theoretical decay amplitude depends on seven complex subtraction constants, one 
of which can be absorbed into overall phase and normalization of the amplitude. The re¬ 
maining parameters are fitted to the experimental Dalitz plot data from the CLEO [21] and 
FOCUS [22] collaborations, restricting the kinematic region to below the r]'K threshold, 
where the elastic approximation is assumed to work well. We have considered different 
ht scenarios with (full) and without crossed-channel rescattering effects (Omnes), as well 
as with and without the ttK isospin 1/2 D-wave. Although the Omnes fits give reason¬ 
able x^/d.o.f., we obtain large destructive interferences between single-variable amplitudes, 
which manifest themselves in unphysical ht fractions. The full hts result in good x^/d.o.f. 
around 1.1 for the CLEO data (1.2 for the FOCUS data), with sensible ht fractions through¬ 
out. Including the ttK isospin 1/2 D-wave does not signihcantly improve the y^/d.o.f., 
however the ht fractions of the nonresonant waves are reduced, giving small interference 
effects between the single-variable amplitudes. We have shown that we can describe the 
^ K~7r~^Tr~^ Dalitz plot data in the region where we deem elastic unitarity to hold 
approximately, solely relying on vrvr and ttK scattering phase shift input and exploiting the 
constraints of dispersion theory. 

Three-body rescattering effects suspends any strict relation between the phase of the de¬ 
cay partial waves and scattering phase shifts: we have shown that the signihcantly stronger 
rise of the ttK 5-wave phases, as observed in analyses of these D-meson decays [23, 26] in 
comparison to phase shift data, can be understood at least qualitatively in the framework 
of Khuri-Treiman equations. 

We have simultaneously constructed the formalism for the decay D~^ —)• K^Tr^TT~^, which 
is directly related to —)• by charge exchange and can be constructed from 

different linear combinations of the same (isospin) amplitudes. This second decay channel 
has recently been measured by the BESIII collaboration [59]. A simultaneous analysis of 
both Dalitz plots will further exploit the predictive power of the dispersive formalism; due 
to the direct contribution of the vrTT P-wave in the Tr'^vr"’' (as opposed to the 7r'’'7r'’') final 
state, we expect to hnd stronger constraints on the subtraction constants featuring directly 
in the corresponding amplitude. The pertinent investigation is in progress [60]. 
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A Inhomogeneities 


In this appendix the inhomogeneities are calculated from Eq. (3.8). To demonstrate the 

7=1 /2 

procedure we will perform the calculation explicitly in the case of ' (s), 

fl^^s) = V3M]/li{s) = K^{s){Tl/\s)+fl/\s)) . (A.l) 


We start with the projection of the decay amplitude Ai _|-+, Eq. (2.5), onto isospin eigen¬ 

states in the s-channel. We introduce the following crossing matrices: 

Mi = Ai/'M/' , Ml = ^ Aii'Mi' (A.2) 

v I' 


and so on, where is the isospin I eigenstate in the x-channel and Xly the crossing 
matrix for the transition from channel y to x, where I and I' are the matrix component 
indices. We obtain the following explicit forms: 



_ 1 / 1 VIo \ 


(A.3) 


The f-channel and rt-channel single-variable amplitudes can be split, with the aid of the 
crossing matrices, into Ig = 1/2 and Ig = 3/2 contributions. 



/s = l/2 /s=3/2 


^l{u) = ^(V3(f - s)J-/(u) + 5J-2 (u)) -liVsit - s)Aliu) - J-2(n)) , 

'-V-''-V-' 

Is=l/2 Is=3/2 

(A.4) 


with L G {0,1}. Retaining the I = lj2 pieces only, we have 




+ [»(* - «) - A] 

+ g (v^(i — s)J-l{u) -|- 5J-'o (u)) . 



(A.5) 


Since there is no isospin 1 component in the ti-channel amplitudes of A4 _|-+, the projections 

onto this specific component yield zero and therefore provide an additional cross-check. 
Similarly no Ig = 1/2 component should appear in A4 qq^. We are left with the angular 
momentum projection. For a compact notation we dehne the angular average integration 
by 

{y)^\ J ^ dzy z^M(x(y, Zy)) . (A.6) 
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We immediately obtain 


, and {z-f). = i-ir {z^f}. 


(A.7) 


The angular average integration is straightforwardly performed in the scattering region. 
The continuation to the decay region, where the naive integration would cross the right- 
hand cut, has been discussed extensively before [6, 15]. We now perform the partial-wave 
projection 

^ j dzsPL{zs)M^^^]!‘^{s,t{s,Zs),u{s,Zs)) , (A.8) 

with the Legendre polynomials Pl{zs). For the S-wave we obtain 

^ + 1 ((-4.2+o.) + 1 [ ( 24 -'" - vioT-y^) 


+ + Bsz + Cs) (2^1 


(A.9) 


where 


_ k(x) 

"" 2x ’ 


a = 


(x^ — A)^ — x^(So — 2x)^ 


Bx = 


Drr = - 


k(x)(x^ -|- A) 

2x2 ’ 

3x2 _ ^ _ 2:5]^ 

2^ 


(A.IO) 


with So = M|, -|- -|- 2M2, X € {s,t}. Thus from Eq. (A.9), the inhomogeneity can be 

immediately read off from the relation = J-’q'^^(s) -|- The full set of 

inhomogeneities is given in terms of the angular averages 






- ( {Auz^ - BuZ - Cu)z^ ( - xl^P't 


2 ^ 3/2 


z[pI'‘^ + 


/i(u) = 


('^o)ns + 2 


Ku{u) 

{AuZ^ - BuZ^ - Cuz) (P^^^ + VTOP^^") 




1 

+ 3 
1 


ipy^ - + ((^2^2 + B,z + C,) (2^^/^ - 


k(s) [ 2 


+|({-4.2^+o.2)x;>. 
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^3/2/ _ 

2 V 6 


1 

~ 6 

2 ^2(s) 


^ -^1 )t 


3 

10 


+ ((^ 2^2 + BsZ + C,) 
oc: c: 

((3.2 - 1) J-2)^^ + _ ((A,.+ 11,) (3.2 - 1) Tl)^^ 

3/2 \ 


2^3 


+ “((32"-i)+''"y -1^ ((322-1)^"; 


° /ts 3 

+ y ((242.2 ^ ^ (3^2 _ jrl / 2 ^ 

where in addition to Eq. (A. 10) we have used 


. 1 / x2 n 1 / \ ^ (^0 — 2n)2 — v? 

Au — ~^^u\^) •) ^ ^ * 


(A.ll) 


(A.12) 


B Invariance gronp matching 


In this appendix, we study the polynomial ambiguities in the decomposition of the total 
decay amplitudes Eq. (2.5) into single-variable functions, dubbed “invariance group”. We 
wish to determine the polynomial at most linear in the Mandelstam variables that can 
be added to the different single-variable amplitudes, leaving the total decay amplitudes 
Eq. (2.5) invariant. For this purpose, we make use of the relation s + t + u = 3so = 
Mq + + 2M^. It is easy to check that adding the following terms to the various 

5-waves as well as the tttt P-wave: 


Tq [u) = ao-\-bou , [u) 

1 /oinv Q /oinv 

(s) = co+do+ rr (^) 


--^bo + 2do 

V5 


2 V 2 


3[oq + 6o(3so ~ 2s)] + 2(co + dgs)^ 


(B.I) 


leaves both Ai _|_+(s,t,u) and Aloo-i-(®)^) unchanged. The most general full decay am¬ 

plitudes are therefore obtained by 


A 


(s) — P^(s) + (s), 


(B.2) 


which, according to Eq. (B.I), has a four-parameter gauge freedom built in. 

Following Ref. [47], we rewrite the polynomial representations of A£^^{s) Eq. (B.I) 
into the Omnes representation ^'£^'{ 3 ) in order to match to Eq. (3.7); 


fmv 


(s) = 0£(s) 



s\n5[{x)F£"'' \ 
\^i{x)\ix -s)i ' 


(B.3) 


with the subtraction polynomials 7r£(s). As the invariance polynomials P'£™'"(s) do not have 
discontinuities, it immediately follows that (s) = which is also conhrmed 
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by a straightforward calculation. We determine the subtraction polynomials by equating 
the polynomial and Omnes representations of the invariance group. We obtain 


7ri(s) = 


T£’^''is) , s" dx sm5i{x)Tj 


ni 


+ 


VT 


.j inv 
L 


\n[{x)\{x-s) 


(B.4) 


The next step is to rewrite the inverse Omnes function into a dispersion relation. Its 
discontinuity is given by 


which thus yields 


disc^ - 

sin (5£(x) 


n' 


1 n POO 

T-^=pUs)-- 


n 

x” |Il£(x)|(x — s) ’ 


(B.5) 


(B.6) 


‘LV 

with the subtraction polynomial PIq{s) = 1 + ■ 4^^® subtraction constants 

(a;£)j are given by the following sum rules, provided that the dispersion integrals converge: 

, r. 1 f°° dx sin (x) 

K)* = -z 

^ Jstt 


Therefore Eq. (B.4) yields 

^i(s) = PLis)pr''{s) + ^ ^ ^ 


^th \^iix)\ 

dx sm6i{x){Tr{x)-Ti^-^s)) 


(B.7) 


\^iix)\{x - s) 


(B.8) 


As an example we will study the single-variable amplitude Pq with = ao + bos. We 

obtain 


'^o(^) ~ + ^0 + oo(wo)i ^ ~ 


1 f°° dx sin5£(x) 


bos 


Using the sum rule value for (cUq)! we find 

= oq + ho + CLo{^o)i 
The other subtraction polynomials are obtained in an analogous way and read 


s. 


(B.9) 


(B.IO) 


7ro(s) = ao + {bo + ao(wo)i)s , 7r((s) = -^^o + 2do , 


-Kq' (s) = co+ do + (a;o^ )i s-F do{ujQ )i-hco(a;o'^ )2 
7^0^^(s) = 36oSo) + 2co 


+ 


do{^y‘^)2 +coiu:^^)^, 


1 / 2 . 


+ 


(a;o^^)i ( V3(ao + 36oSo) + 2co ) - 2 ( V35o - 


(B.ll) 


with no contributions to the ttK P- and D-waves. Polynomial terms with higher order 
than the subtraction polynomials of the corresponding amplitudes (see Sec. 3.3) have been 
omitted. 
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As we have argued above that a choice of the constants qq, ■ ■ ■, dg corresponds to a 
mere “gauge” choice and is unobservable, we can decide to fix them by requiring the (linear) 
subtraction polynomials in the nonresonant S'-waves (1 = 2 tttt and 1 = 3/2 ttK) to vanish. 
Equation (B.ll) proves that this is feasible: we can eliminate all subtraction constants in 
J'q by the appropriate choice of oq and bg, and all constants in by adjusting cg and 
dg. The result is the system Eq. (3.10) in the main text, which is thus free of ambiguities. 

C Numerical treatment 

In this appendix we discuss the numerical treatment of the double integral in Eq. (3.13), 
i.e. the part independent of the subtraction constants. In the following we will restrict 
ourselves to the s-channel case for illustration. We rewrite the term in Eq. (3.13) as 
z^(s, t) = C"^{s, t)/K^{s), with (/(s, t) = (2ts — 3sos -|- — A) and 3so = -|- 2M^, 

such that the double integral adopts the form 




S K 


2L—m 



F(x) sin(5(x) 
x'^\U.(x)\k?‘^^^ (x) 



X — t 


dt dx, 


(C.l) 


with t±(s) = t(s,±l). Eirst we study the case s > (Md — Mt^)^. The angular integral can 
directly be performed as the two integral paths do not cross each other. We may simply 
use 





T{x) sin(5(x) 
X"’|n(x)|K2'^+l(x) 


W(s,x) dx, 


W{s,x) = 


(s) ^ ^ 


dt , 
(C.2) 


where W(s,x) can be determined numerically in a straightforward way. The discretized 
integral reads 



T(x) sin 6 (x) 
x"-|n(x)|«;2^+^(x) 


IE(s,x)dx 



(C.3) 


where Cq(s) -t-c((s)x is the linear interpolation of VE(s, x) sin (5(x)/x”|fl(x)| in the interval 
[sj,Sj-|_i] for a fixed s. Note that the resulting integrals can be performed analytically 
with the singularities moved into the upper complex plane to obtain the correct (physical) 
branch. 

Eor the case s < (Md — M^^)^ the Cauchy singularity needs to be handled carefully, as 
the integration paths meet. We rewrite 


poo 

/ 

j S+Vi 


sin(I(x) 


/Sth X^\n{x)\K 

/"(x) sin(I(x) 


2L+1 


L 


(x) Jt-is) 




T(x)sm6{x) r 

L x-|I2(x)|k2^+i(x) U ^ ^ 


^^dx 

t_ (s) ^ ^ 

rn(t)-x^n{x) 


(s) 


x — t 


dt dx 


^ .oo ri+(^) 

J Sth 


t_{s) K(x)2i+l(x-t) 


(C.4) 
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The first summand is treated as above in Eq. (C.3). For the second summand we obtain 


OO ft+is) 


poo P 

Jsth Jt- 


sin(5(x) 

't-{s) -t) 


dt dx = 




ET»y) rv 

J Sn •J t. 


Oq + a\x 


t-(s) l^'^^+^{x){x-t) 


dt dx, 


(C.5) 


where now Og + a\x is the linear interpolation of e**^^*^ sin(5(x) in the interval [sj, Sj+i]. 
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